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ABSTRACT 

We carried out the direct measurement of the fraction of dusty sources in a 
sample of extremely red galaxies with (R — K s ) > 5.3mag and K s < 20.2mag, 



> 

<^j- ; using 24/im data from the Spitzer Space Telescope. Combining deep 24/im, K, 



and i?-band data over an area of ~64 sq.arcmin in ELAIS Nl of the Spitzer First 



m 
\o 

Look Survey (FLS), we find that 50±6% of our ERO sample have measurable 
24/im flux above the 3cr flux limit of 40/iJy. This flux limit corresponds to a SFR 
of 12M /yr at z ~ 1, much more sensitive than any previous long wavelength 
measurement. The 24/im-detected EROs have 24-to-2.2 and 24-to-0.7/im flux 
ratios consistent with infrared luminous, dusty sources at z > 1, and an order 
of magnitude too red to be explained by an infrared quiescent spiral or a pure 
old stellar population at any redshift. Some of these 24/im-detected EROs could 
be AGN, however, the fraction among the whole ERO sample is probably small, 
10-20%, as suggested by deep X-ray observations as well as optical spectroscopy. 
Keck optical spectroscopy of a sample of similarly selected EROs in the FLS field 
suggests that most of the EROs in ELAIS Nl are probably at z ~ 1. The mean 
24/im flux (167/iJy) of the 24/im-detected ERO sample roughly corresponds to 
the rest-frame 12/im luminosity {vL v {12^im)) of 3 x 1O 1O L at z ~ 1. Using the 
correlation between IRAS vL v {Yl[im) and infrared luminosity L/_r(8 — 1000/i), 
we infer that the < L IR > of the 24/im-detected EROs is 3 x lO n L and 1O 12 L 
at z — 1.0, 1.5 respectively, similar to that of local LIRGs and ULIGs. The 
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corresponding SFR would be roughly 50 — 17OM /yr. If the time scale of this 
starbursting phase is on the order of 10 8 yr as inferred for the local LIRGs and 
ULIGs, the lower limit on the masses of these 24/zm-detected EROs is 5 x 10 9 — 
2 x 1O 1O M . It is plausible that some of the starburst EROs are in the midst 
of violent transformation to become massive early type galaxies at the epoch of 
2-1-2. 

Subject headings: galaxies: bulges - galaxies: starbursts - galaxies: infrared 
luminous - galaxies: galaxy evolution - galaxies: high-redshifts 

1. Introduction 

Optical/near-IR colors, such as (R — K s ) or (I — K s ) have been commonly used in wide- 
area surveys to select old stellar populations at z ~ 1 — 2 (Cimatti et al. 2002; McCarthy et 
al. 2001; McCarthy 2004 for the review of this subject). Near-IR observations, which sample 
cool low-mass stars, are sensitive to old stellar populations. The (R — K s ) of 5.3 corresponds 
to the calculated color of a passively evolving elliptical galaxy at z — 1. Therefore, in 
principle, the color criterion of (R — K s ) > 5.3mag or (7 — K s ) > 4mag (EROs) should 
select early type galaxies at 2 ~ 1. However, these color selections are also sensitive to 
dust-reddened, star-forming systems, and examples of both passively evolving ellipticals and 
dusty starburst EROs have been found (Soifer et al. 1999; Hu & Ridgway 1994; McCarthy, 
Persson & West 1992; Graham & Dey 1996). This indicates that the optical/near-IR SEDs of 
these sources are sufficiently degenerate that these color criteria cannot effectively distinguish 
between them. In addition, a small fraction of EROs (10-20%) could also be AGN, as shown 
by deep Chandra data and optical spectroscopy (Alexander et al. 2002; Yan, Thompson & 
Soifer 2004). The relative contribution of these two galaxy types — old stellar populations 
and dust-reddened, star-forming galaxies — is a critical issue for many surveys whose goal 
is to determine the evolution of the mass function and the formation of massive galaxies. 
Deep optical spectroscopy indicates that a large fraction (30-50%) of EROs have emission 
lines (Cimatti et al. 2002; Yan, Thompson & Soifer 2004; McCarthy et al. 2004); however, 
it remains unclear what fraction of EROs are truly dust obscured galaxies. 

MIPS 24/im data from the Spitzer Space Telescope (Rieke et al. 2004; Werner et al. 
2004) offer the first opportunity to directly address this critical issue. Dusty, star-forming 
galaxies are clearly distinguished from early-type galaxies at mid-IR wavelengths. Between 
z ~ 1 — 2, the MIPS data is especially discriminating as strong, rest-frame 6— 12yU.ni polycyclic 
aromatic hydrocarbon (PAH) dust features redshift into the 24/im-band. In this Letter, we 
present our initial study of the 24/im properties of the (R — K s > 5.3)mag selected EROs in 
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ELAIS Nl. Throughout the paper, we adopt H = 70km/s/Mpc, tt M = 0.3, f2 A = 0.7, and 
the Vega system for optical/NIR magnitudes. 

2. Data 

2.1. Spitzer 24/im Observations and Data Reduction 

The primary dataset used in this Letter is in the ELAIS Nl field, which is a part of the 
Spitzer First Look Survey (FLS) 1 . The field was observed in a 2 x2 photometry mode mosaic, 
and the raw data were processed and stacked by the data processing pipeline at the Spitzer 
Science Center (SSC). Source catalogs at 24/im were generated using StarFinder (Diolaiti 
et al. 2000), which measures profile-fitted fluxes for point sources. The complete description 
of the 24/xm data reduction and source catalog can be found in Fadda et al. (2004b) and 
Marleau et al. (2004). To aid in the interpretation of the results from ELAIS Nl, we also 
analyzed the FLS verification strip (FLSV). The data presented here covers 64 sq.arcmin in 
ELAIS Nl and 256 sq.arcmin in the FLSV. Table 1 summarizes the salient characteristics of 
all data used in this paper. 

2.2. Optical, Near-IR Imaging and Keck Spectroscopy 

All .R-band observations were taken at the Kitt Peak National Observatory. Final 
stacked images and source catalogs have been publicly released (Fadda et al. 2004a). The 
fC,-band data were obtained using the Wide-Field Infrared Camera (WIRC) on the Palomar 
200-inch telescope. The K — band data in Elais Nl covers 8' x 8'. A detailed description 
of the WIRC observations and data reduction will included in a separate paper by Choi et 
al. (2004). Finally, high- resolution, optical spectra of a ^-selected galaxy sample within 
the FLS region were obtained using the Keck, Deep Imaging Multi-Object Spectrograph 
(DEIMOS) (Faber et al. 2003). A total of « 1000 redshifts were measured, of which 112 (52 
EROs) are included in the analysis in this paper. 



For details of the FLS observation plan and the data release, see http://ssc.spitzer.caltech.edu/fls. 
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3. Results and Implications 

3.1. The 24/xm detected EROs 

To merge the R/K and 24/im source catalogs, we used a simple positional matching 
method with a 2.4" match radius, which corresponds to 3a combined astrometric uncertainty 
from the K s and 24/xm data. Due to the relatively low source density (7arcmin~ 2 ), the 
likelihood of spurious matches is small (3%), consistent with the fact that we find only one 
multiple match out of 65. In addition, we test the robustness of each match by comparing 
the probability of the measured separation based on the astrometric uncertainties to the 
probability of a spurious detection based on the 24/im source density. The probability ratios 
above unity imply to be likely real matches, and we found that all but two matches meet this 
criteria. We retained these two matches since the visual inspection suggests that they could 
still be associated with physical offset centroids. Bright stars were rejected using i?-band 
images. Stellar contamination in our sample is expected to be small since the field is at 
the galactic latitude of 41°, and 24/xm data samples the tail of the Rayleigh- Jeans energy 
distribution. 

Since the R- and JT s -band data have similar seeing, (R — K s ) colors were measured using 
a fixed 3" diameter aperture. The ERO catalog with (R — K s ) > 5.3mag and K s < 20.2mag 
(6a) consists of 129 galaxies over 64 sq.arcmin in ELAIS Nl. Figure 1 shows the (R — K s ) 
vs. K s distribution for all sources (black crosses. The (R — K s ) > 5.3mag and K s < 20.2mag 
limits are shown as solid horizontal and solid vertical lines, and sources with 24/im detected 
counterparts are indicated as red open circles. We find that 24/im detected sources have 
slightly redder (R — K s ) colors than non-detected sources, consistent with the expectation 
that 24/im emission is an indicator of dust extinction. 

Of the 129 EROs, 65 (50±6%) have 24/xm emission with flux greater than 40/iJy. The 
fraction of 24/xm detected EROs becomes slightly higher for redder sources, with 56±15% 
and 67±30% for (R — K s ) > 6.0mag and (R — K s ) > 6.5mag respectively. However, the 
errors of these fractions are large due to the shallow R band limit. Deeper data would be 
needed to reduce the uncertainties. In figure 2a, we show the 24/xm flux distributions of both 
the total and ERO samples. We find that 70% of the 24/xm EROs have iW^m > 90/xJy and 
that the mean 24/xm flux of the ERO sample is 167/iJy. Figure 2b shows K s vs. 24/xm flux 
for all detected sources (solid squares) and for the EROs (large open circles). We find that 
amongst the 24/xm detected sources, K s is weakly correlated with /24/im, with large scatter. 
This is expected since these two bands sample light from different physical origins, stellar 
photosphere versus dust emission. This also explains why many faint K s sources have fairly 
bright 24/im fluxes. For K s bright (~ 18 — 17) and 24/im-faint sources, the mid-infrared 



data is deep enough to detect the emission from normal galaxy populations at low redshifts. 
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Fig. 1. — The (R — K s ) vs. K s color- magnitude diagram. All i^<,-band detected sources are 
shown as crosses, and sources with 24/im counterparts are indicated by open circles. The 
(R — K s ) > 5.3mag and K s < 20.2mag limits are shown as solid horizontal and solid vertical 
lines. 
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Fig. 2. — Panel 2a shows the 24/zm flux distribution of the total (dashed line) and 24/im- 
detected ERO (solid line) populations. Panel 2b shows K s vs. 24/im flux for all detected 
sources (solid squares) and for the EROs (large open circles). 



3.2. The Nature of the EROs with 24/im Emission 



What types of galaxies are the 24/zm detected EROs? In Figure 3a, we present a 
24-to-0.7//m and 24-to-2.2/xm color-color diagram showing that 24/zm-detected EROs (red 
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squares) are clearly separated from EROs without 24/im emission {green points), as well 
as the general non-ERO population {black and blue). This suggests that the MIR/Optical 
and MIR/NIR colors of 24/im-detected EROs are unique, distinguishing them from other 
populations. Figure 3b presents the expected colors as a function of redshift for various SED 
templates, assuming no evolution. The SEDs for M51 (normal spiral), M82 (starburst) and 
Mrk231 (dusty AGN) are taken from Silva et al. (1998), Fadda et al. (2002), and Chary & 
Elbaz (2001). In the case of M31's bulge (old stellar population), we use the near-IR J, H, 
K, IRAS 12, 25, 60 and 100/xm fluxes, all within a 4' diameter aperture of the central nucleus 
(Soifer et al. 1986) and merge this with a 10 Gyr theoretical SED from Bruzual & Chariot. 
2 For comparison, a sample of 112, 24/im sources from the FLSV with known spectroscopic 
redshifts are marked to show the actual redshift range of the 24/im EROs in ELAIS Nl. 
Of these 112 redshifts, 10 are 24/im-detected EROs with 0.8 < z spec < 1.3, and the optical 
spectra of these 10 sources all have [OII]3727A emission line. 

Comparing Figure 3a & 3b, we reach the following conclusions: 1) 24/im-detected EROs 
in ELAIS Nl with (R - K s ) > 5.3mag, K s < 20.2mag, and /„(24//m)> 40/iJy are infrared 
bright sources at z > 1. They have colors similar to starbursts like M82 at z > 1, or dust 
reddened AGN like Mrk231 at z > 0.7. Their colors are too red to be explained by any 
normal spiral or old stellar populations at any redshifts. 2) The likely redshift range of 
these 24/im EROs is z > 1, as predicted from the model SEDs of M82 and Mrk231. This is 
further confirmed by comparison to the Keck spectroscopic sample from the FLSV. 3). The 
remaining half of the ERO population are probably galaxies with old stellar populations at 
z ~ 1, as suggested by the tracks in Figure 3b. 

Half of the ELAIS Nl ERO sample have 24/im dust emission. The key question is 
what their infrared luminosities are. In the previous section, we conclude that the 24/im- 
detected EROs have colors similar to M82, however, Figure 3a does not set any constraints 
on either their luminosities or masses. We have a total of 112 redshifts for objects with 
24/im counterparts (the FLSV 3a flux limit is 90/iJy). In Figure 4, we present the observed 
24/xm luminosity versus redshift for the spectroscopic FLSV sample. At z ~ 1 the observed 
24/im luminosity roughly corresponds to the rest-frame IRAS 12/mi luminosity, within a 
factor of 2. The Spitzer 24/im filter is narrower than the IRAS 12/im filter, but ignoring this 
difference, a crude conversion of the observed Spitzer 24/im flux can be made to place a lower 
limit on the rest-frame IRAS 12/im luminosity. We can infer the total infrared luminosity of 
24/im sources at z ~ 1 by using the correlation between the 12/im luminosity {yL v {Yl\xm)) 
and the infrared luminosity L IR (8 - 1000/im), L IR = 0.89 x (u L v (12 /im)) LO94 L (Soifer 



2 ftp: / / gemini.tuc.noao.edu/pub/charlot /bca5 
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et al. 1989; Chary & Elbaz 2001). The mean flux (167// Jy) of the 24/tm-detected EROs 
in ELAIS Nl implies the L IR ~ 3 x 10 n L o at z = 1.0 and L IR ~ 1O 12 L at z = 1.5, 
similar to that of local LIRGs and ULIGs. The corresponding SFR is 50 — 170M Q /yr, using 
SFR = 1.71 x 10- 10 L IR (M Q /yr) (Kennicutt 1998). We emphasize that the 3a flux limit 
of 40/iJy in ELAIS Nl corresponds to a SFR of 12M /yr, more sensitive than the deepest 
1.4GHz observation (Smail et al. 2002). In comparison, the SFR derived from the rest- 
frame [OII]3727A emission line for EROs is roughly a few xM /yr (Yan, Thompson & Soifer 
2004; Cimatti et al. 2003). Our results are consistent with very deep 1.4GHz (rms ~ 3/tJy) 
measurements (Smail et al. 2002). 



3.3. Summary & Implications 

Spitzer 24/tm data in ELAIS Nl have revealed that 50±6% of EROs with (R - K s ) > 
5.3mag and K s < 20.2mag have detectable 24/tm emission above 40/tJy. The colors and 
inferred redshifts of the 24/tm-detected EROs suggest that they are infrared- luminous, dusty 
sources at z > 1. Their mean 24/tm flux (167/tJy) corresponds to < L IR >~ 3 x 10 n L Q 
at z ~ 1 and to < L IR >~ 1O 12 L at z ~ 1.5. Some of these 24/tm-detected EROs could 
be AGN. The one mega-second Chandra observation in the HDF and optical spectroscopy 
suggest that the fraction of EROs likely to be AGN is small, around 10-20% (Alexander 
et al. 2002; Yan, Thompson & Soifer 2004). These dusty AGN can be identified when we 
combine this current analysis with the IRAC data. Our result suggests that a significant 
fraction of EROs are extremely active starbursts (LIRGs or ULIGs). If the time scale of 
this starbursting phase is on the order of 10 8 /yr, as inferred from local LIRGs and ULIGs 
(Sanders & Mirabel 1996), we can set a lower limit to the mass of these 24/tm-detected EROs 
as SFRxr = 50 - 170 x 10 s = 5 x 10 9 - 2 x 1O 1O M . If the EROs without detectable 24/tm 
are indeed massive systems with old stellar populations at z ~ 1 as measured by several 
recent surveys (Glazebrook et al. 2004; Bell et al. 2003), one plausible connection between 
the starburst and early-type ERO populations is that the former may be in the process of 
transforming into the latter, as initially postulated by Kormendy & Sanders (1992). In the 
hierarchical clustering paradigm, it could be interpreted that our deep 24/tm observations 
are capturing a massive galaxy population in the midst of violent transformation - possibly 
in the process of assembly via mergers/starbursts at the epoch of z ~ 1 — 2. The accurate 
determination of the stellar and dynamical masses of these starburst EROs at z ~ 1 — 2 will 
be critical for the resolution of this question. Finally, the measurement of volume averaged 
mass density at z ~ 1 — 2 would require a better understanding of the physical source of the 
integrated K-band light — whether from dusty systems or from old stars. 
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Our result is consistent with what has been found with HST morphological studies of 
EROs. Yan & Thompson (2003) and Moustakas et al. (2004) have found that close to 50% 
of EROs with K s < 20 have morphologies consistent with disk or later type galaxies in 
the observed 8100-8500 A wavelength, and less than 40% show clean bulge type profiles. 
With the HST/ACS/NICMOS images in the FLSV region, we will be able to investigate 
the morphologies of these infrared luminous EROs, and to determine if indeed they are 
starbursting mergers at z ~ 1. 

This work is based in part on observations made with the Spitzer Space Telescope, 
which is operated by the Jet Propulsion Laboratory, California Institute of Technology under 
NASA contract 1407. Support for this work was provided by NASA. The spectroscopic 
data presented herein were obtained at the W.M. Keck Observatory, which is operated 
as a scientific partnership among the California Institute of Technology, the University of 
California, and the National Aeronautics and Space Administration. The Observatory was 
made possible by the generous financial support of the W.M. Keck Foundation. We also 
wish to recognize and acknowledge the very significant cultural role and reverence that the 
summit of Mauna Kea has always had within the indigenous Hawaiian community. We are 
most fortunate to have the opportunity to conduct observations from this mountain. 
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Fig. 3.— The two panels show log 10 (f/„(24//.m) / uf u (0.7/im)) versus 

log 10 (uf v (2Afim)/uf v (8fim)). Panel 3a presents the data from ELAIS Nl. In this 
panel, red/ 'green squares represent EROs with/without 24/zm detections. For the non-24/im 
EROs, the green crosses have only upper limits in the i?-band. The black/blue points 
represent the full ^-selected sample with/without 24/im detections. The plotted colors of 
sources not detected at 24/xm are upper limits and could move down to lower-right corner of 
the diagram, as illustrated by the green & blue arrows. The black solid line in both panels 
represents constant (R — K s ) = 5.3mag color. Panel 3b shows the computed color-color 
tracks as a function of redshift for a set of SED templates of local, well-known galaxies. 
The subset of the FLSV sample with measured spectroscopic redshifts are overplotted for 
three different redshift bins: z < 0.5 (magenta open triangles); 0.5 < z < 1.0 (orange open 
squares); and z > 1.0 (cyan open circles). 
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Fig. 4. — The observed 24/im luminosity as a function of spectroscopic redshift for a sample 
of 112 24/zm-detected sources in the FLSV region. The solid triangles are 24/zm-detected 
EROs. 
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Table 1. Summary of Optical, NIR and 24/im Observations 



Filter 


R 


K s 


24/^m 


FWHM (//) 


1.0 


1.0 


5.5 


Pixel Scale (///pixel) 


0.26 


0.25 


2.55 


<Exp. Time> (sec) 


1800 


1800 a /9000 b 


349 a /4268 b 


3<r Limits 


25. 5 a /25.2 b 


20.3 a /21.0 b 


90 a /40 b 




mag 


mag 





a FLSV 
b ELAISNl 
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